Josephson-current induced conformational switching of a molecular quantum dot 
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We discuss the behavior of a two-level system coupled to a quantum dot contacted by super- 
conducting source/drain electrodes, representing a simple model for the conformational degree of 
freedom of a molecular dot or a break junction. The Josephson current is shown to induce confor- 
mational changes, including a complete reversal. For small bias voltage, periodic conformational 
motions induced by Landau-Zener transitions between Andreev states are predicted. 

PACS numbers: 74.50.-|-r, 74.78.Na, 73.63.-b 



The remarkable recent progress in the fabrication 
and experimental study of transport through ultrasmall 
nanoscopic devices, break junctions, or molecules (in the 
following termed "quantum dot" or simply "dot" ) has 
stimulated renewed interest in the Josephson effect [J, 
where the Josephson current through a dot contacted by 
superconducting electrodes with phase difference ip is the 
relevant observable. The full current-phase relation has 
been measured in various systems, and electron-electron 
interactions on the dot were shown to be important 
as expected from theory j^]. The already achieved wide 
tunability (via gate electrodes) and impressive control 
over Josephson currents through nanoscalc dots indicate 
that experiments should be able to also probe modifi- 
cations of the super-current due to the coupling of the 
dot to another quantum system (e.g., a spin or a side- 
coupled dot). Many previous efforts have focussed on 
studying the coupling to the spin degree of freedom in 
molecular magnets Q, which is also related to issues ap- 
pearing for superconductor-ferromagnet-superconductor 
structures . Theoretical work has also discussed the ef- 
fects of local vibration modes on the super-current .where 
the dot is coupled to a boson mode (phonon) 0, Q . 

Surprisingly, so far the effects of a two-level system 
(TLS) coupled to the dot have not been addressed, ex- 
cept for normal-conducting leads. This is an important 
question, since for instance two conformational conGgura- 
tions of a molecule may represent the TLS degree of free- 
dom. Experimental results for molecular dots or break 
junctions (with normal leads) were interpreted along this 
line [1, [l^, [ll|, [l^ , but a TLS can also be realized for a 
side-coupled double-dot system in the Coulomb block- 
ade regime For concreteness, we here refer to the 
TLS states CTz = ±1 as the two distinct conformational 
states of a molecular dot, where Cz couples to the dot's 
charge. A coupling of the TLS to the dot's spin does 
not have a significant effect on the phenomena of inter- 
est here, see below. (In any case, spin effects have been 
addressed in different contexts previously Q.) Our the- 
ory indicates that by variation of the phase if, the TLS 
state can be significantly affected over a wide parameter 
regime, including a complete reversal of the conforma- 
tional configuration. This remarkable effect allows for 
the dissipationless control (including switching) of the 
conformational degree of freedom (cr^) in terms of the 



phase difference ip, which can be tuned experimentally 
by embedding the device in a SQUID geometry Q. Con- 
versely, changing the conformational state will affect the 
Josephson current in a distinct manner. Moreover, when 
applying a bias voltage, a periodic conformational mo- 
tion is triggered via the ac Josephson effect involving 
Landau-Zener (LZ) transitions between Andreev states. 
Our predictions (both for zero and finite bias) can be 
tested experimentally for a wide class of molecules elec- 
trically contacted in a break junction setup. Related ex- 
periments, reporting TLS behavior due to a conforma- 
tional variable, have been published for normal leads 0]. 
For normal leads, the model employed below has also 

Avail- 



been motivated in a recent theoretical work 
able parameter estimates for dot and TLS energy scales 
suggest that the predicted phenomena can be ob- 
served using existing state-of-the-art experiments. De- 
tection schemes to read out the conformational state are 
also available, e.g. by single-molecule force microscopy 

H. 

We study a spin-degenerate molecular dot level with 
single-particle energy ed and on-site Coulomb repulsion 
U > 0, coupled to the TLS and to two superconduct- 
ing banks (leads). Employing the standard wide-band 
approximation for the leads, we assume a symmetric sit- 
uation [l^, where the banks arc modelled as identical 
s-wave BCS superconductors with gap A and the dot- 
lead hybridizations are equal, = Tji = T/2. The TLS 
describing the conformational state corresponds to Pauli 
matrices cr^.z, with bare energy difference Eq and tunnel 
matrix element Wq between the two states. The Hamil- 
tonian is H = Hq + iJtun + ^^loads, where the coupled 
dot-plus-TLS part is (we set e ^ h ^ ks = I) 
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: dld^ for the dot 
describes standard 



with the occupation number Ug = 
fermion dg with spin s =t,i, -fficads 
BCS Hamiltonians, tp can be included by phase factors in 
the tunnel Hamiltonian fftun , and we define the renor- 
mahzed dot level e = ed+U/2. In Eq. ((T|), the TLS couples 
with strength A to the charge on the dot, which can be ra- 
tionalized in simple terms by assuming a one-dimensional 
effective reaction coordinate X describing conformations 
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of the molecule. The dominant coupling to the electronic 
degrees of freedom is then (as for phonons) of the form 
(X + ni) 0, S, 0, In the Hmit of interest, the 

potential energy V{X) is bistable with two local minima, 
and a truncation of the low-energy dynamics of X to the 
lowest quantum state in each well leads to Eq. ([T]). For 
a detailed derivation, see also Ref. [l^ . 

When dealing with the equilibrium problem, it is con- 
venient to work with Nambu spinors d^r) = [d-^ , d| ) in 
imaginary time r. The lead fermions can then be inte- 
grated out exactly, and the partition function is Z = 

Tr (^e--f^"/^re-/''^''^'''*(^)^(^~^')''(^')^ , where the trace 

extends over the dot-plus-TLS degrees of freedom only, 
T denotes time-ordering, and the effect of the BCS leads 
is contained in the 2x2 Nambu self-energy matrix E(t), 
whose Fourier transform is |4|1 



E([j;^) 
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Acos(i^/2) \ 
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We mostly consider zero temperature, T = 0, where both 
the Josephson current /(</?) through the dot and the ex- 
pectation value S = (cTj,) of the conformational state fol- 
low from the ground-state energy Eg{ip^ Eq) according to 



2^. 
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(3) 



Later on the formalism will be extended along the lines 
of Refs. [IB, to allow for the description of a small 
bias voltage V as well. 

Let us first illustrate our central findings when both 
the charging energy U and the tunnel splitting Wq are 
very small. Later on we show that for sufficiently small 
U <Uc, finite U has no effect. The ground-state energy 
Eg = \Tm\{E+^E-) then follows from the energies E^r = 
a{X — Eq)/2 — e-^{ip) for fixed conformational state tr = ± 
with dot level = e + o-A/2. With Eq. Q and = 
diag(l, —1), the Andrccv state energy for arbitrary A/F 
follows from 



J 27r det [uv — r^e^ — 2^[uj; U)\ 
In the hmits F > A and A > F, this yields [2] 



= A, 



l-T,sin2(^/2), 
^ F » A, 



(4) 
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with the normal transmission probability 7^ = [1 + 
e^/F^]-!. As long as E+ < E^ [E^ < E+), we have 
S{lp) = H-l(— 1), i.e. the conformational state tr = ) 
is realized, with ideal (perfect) switching when the bands 
E^{ip) and E^{ip) cross at some phase < ip* < tt. 
Hence a necessary condition for switching follows: one of 
the two inequality chains (with TZ^ = 1 — 7^ ) 



must be obeyed. If the dot level is close to a resonance, 
e+ « or e_ « 0, the reflection probabihties and TZ- 
arc significantly different, and Eq. ([6]) holds over a wide 
parameter range. Then Eq. ^ yields 
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In the regime (jB]), the transmission amplitude switches 
between and 7I_ when cp ~ ip* . This implies non- 
standard current-phase relations, as shown in the upper 
inset of Fig. [TJ 

Having established the basic phenomenon, we now 
address the effects of finite U and/or tunneling Wq. 
Progress can be made in the limits F A and A ^ F. 
Let us start with the case when A is the largest energy 
scale of relevance. Then the dynamics is always con- 
fined to the subgap regime (Andrecv states), and quasi- 
particle tunneling processes from the leads (continuum 
states) are negligible. Technically, Eq. ([2]) can then be re- 
placed by S(t) = r cos{p/2)S{t) ^ ^ ^ ^ , and the prob- 
lem is equivalently described by the effective Hamiltonian 
Hcff = i?o -I- Fcos(.^/2) + djdjy The resulting 

Hilbert space can be decomposed into orthogonal sub- 
spaces, TC = TCa (8) TCs, where the Andrecv sector Ha is 
spanned by the zero- and two-clcctron dot states |0) and 
|2) ~ djdjlO) (and, of course, by the conformational TLS 
states), while Tis is spanned by the one-electron states 
|s) = dl\0). For convenience shifting H^s H^g — e, 
the single-particle sector has a pair of d oubly-degenerate 
eigenenergies — ± eg with es = ^\/{Ea — A)^ + Wq, 
whereas the Andreev sector is described by 



\ — En 



Wq 
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Eq^A^ 



A_ 



(6) 



+Tcoa{p/2)T^ 
(8) 

with Pauli matrices t^^z acting in {|2),|0)} subspace. 
If the groimd state of Hcs lies in the Andreev sector, 
the Josephson current can be non-zero, while otherwise 
7 = due to the i^s-indepcndence of the single-particle 
sector. For sufficiently strong interactions, U > Udp), 
the ground state of HcS is in the single-particle sector 
Ti.s- This is indicative of a quantum phase transition to 
the magnetic 7r-junction regime [1]. While this regime is 
outside the scope of Eq. ([8]) (since continuum states are 
not included), we have confirmed this scenario by a per- 
turbative calculation expa nding in F for the full model. 
For A ^ 0, we find Uc = 2^e2 + F2 cos2(.^/2), see Eq. © 
for A ^ F. Note that e and hence Uc can in principle 
be tuned by a gate voltage. For A » maxdi^ol, |e|,r), 
we instead find Uc = A. Because is independent 
of U (up to the shift e = + U/2), a weak interaction 
U < Uc has no effect, and in what follows we set U = 0. 
Since a coupling of the TLS to the dot's spin involves 
only the (^-independent subspace Tis, such couplings are 
also of little relevance for switching, in accordance with 
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the small polar displacements predicted for spins in a 
Josephson junction 

Physical observables can then be computed from 
in Eq. ([5]). The eigenenergies are roots to the exactly 
solvable quartic equation 



E'' - 2A,S^ + AiS + A 







0, 



(9) 



with coefficients A2 = + e| + Ai = 2\e{Eo - A) 

and Ao = {e\ - e| + \^ /^f - X^{e^ /i). The lowest- 
lying of the four roots yields the exact but lengthy result 
for the ground-state energy Eg. Convenient expressions 
for S{(p) and I{ip) in Eq. ([3]) follow by taking the respec- 
tive derivatives directly in Eq. ([9]). For instance, with 
A^ = dAi/dEo, the conformational variable reads 



2K'^El 



A[Eg - A', 



2Eg{El-A2) + Ai/2' 



(10) 



Typical results for S{lp) and are shown in Fig. [T] 
The most efficient way to induce conformational changes, 
including a complete (symmetric) reversal S — S", is 
achieved in the weak-coupling regime A ^ e^, where the 
four roots to Eq. Q can be simplified to 



(11) 

with ground-state energy Eg = E Remarkably, 

Eq. (jlip remains accurate even for A « e^i. A com- 
plete reversal is achieved when tuning Eq or e such that 
E„ = A[l-sgn(e)]-.F with.^ = -|sgn (e) [1 - |£| /eA(0)]. 
In that case, S{0) = -S{tt) = J^/y/wfT^. When 
comparing to the Wq — result, we observe that a fi- 
nite tunnel amplitude Wq only leads to a rounding of the 
transition and a decrease in the switching amplitude, but 
it does not destroy the effect. Finally, with Eq. (fTTj) . the 
Josephson current in the weak-coupling limit is 



sin (f 



Xe 



S{cp) 



(12) 



Next we briefiy discuss the opposite limit within a sim- 
ilar truncation scheme, setting U = 0. For F ^ A and 
if ^ 27rn (integer n), the relevant subgap dynamics is 
again captured by an effective two-level Hamiltonian de- 
scribing the Andreev states [13 , coupled to the conforma- 
tional TLS. With Pauli matrices r^.y^z in Andreev level 
subspace and the notation [see Eq. ([5])] 

= A,e-"y^^/^ sin((^/2)r, + cos(y./2)7 

the effective Hamiltonian follows in TLS space as 



cff 



A-_Bn 



Wo 
2 



Wo 
2 



(13) 



Physical observables are then easily obtained, see the 
lower inset of Fig. [TJ Again the qualitative features of 
the Wq = solution persist. 
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FIG. 1: Conformational state S{(p) — S'(— vs super- 
conductor's phase difference ip. Results from Eq. (|10|) for 
A ^ r are shown for tunnel amplitudes Wo = (dotted) and 
Wo = 0.04r (solid), with A = e = r/2 and Eo = 0.14r. The 
dashed curve gives the exact result for Wq = and A = SF, 
see Eq. Q, extended to finite temperature T = O.OIF. The 
upper inset shows the corresponding Josephson current-phase 
relations. Lower inset: Same as main figure but for F — 4A 
with A = 2e = A/2 and Eo = 0.45A. The dotted (sohd) curve 
is obtained from the F ^ A effective Hamiltonian (|13|l with 
VFo = (VFo = 0.04A). The exact resuff for VFo = is shown 
as dashed curve for T = O.OIA. 




FIG. 2: Time dependence of S (upper row) and ac Josephson 
current / (lower row). The left panel shows the adiabatic 
evolution for low voltage, V — O.OIF. Parameters are the 
same as for the solid curve in the main panel of Fig. [1] The 
right panel is for 1/ = 5F with e/F = 0.2 (sohd) and 0.6 
(dashed, current not shown). Other parameters are A = r/2 
and Eo = Wo = 0.2F. 



The effective Hamiltonian ([S]) for A ^ F also allows 
to study the voltage-biased junction with V" ^ A, where 
the superconducting phase difference is time-dependent, 
ip{t) = 2Vt. During the time evolution induced by ip{t), 
the Andreev and single-particle Hilbert subspaces Ha 
and Hs remain decoupled and mutually orthogonal. The 
task is therefore reduced to solving the time-dependent 
Schrodinger equation i9t*(t) = iJ^(t)*(t), where *(t) 
is a 4-componcnt wave function representing the two An- 
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drccv states and the TLS, and H^{t) is given by Eq. ^ 
with (f 2Vt. For this description to hold at finite 
A, the escape rate 7 of Andreev state quasiparticles into 
the continuum states of the leads should be negligibly 
small. The rate 7 follows from the tunneling self-energy, 
see Ref. [13] for the opposite limit F S> A. For e = 0, we 
find 

7^rexp(^-^[ln(2A/F)-l]), (14) 

leading to exponentially small rates for realistic system 
parameters throughout the regime A ^ F. Numerical so- 
lution of the time-dependent Schrodinger equation leads 
to the results in Fig. [21 They can be understood in terms 
of the four eigenenergies p4|) . For small V, the time evo- 
lution is basically adiabatic, and the LZ probability is 
very small. The left panel in Fig. [2] shows such an adia- 
batic evolution involving time-periodic level crossings of 
the bands and E |_ in Eq. (fTTI) . thereby explain- 

ing the existence of two different supercurrent oscillation 
amplitudes. The "noisy" features in S{t) are fully repro- 
ducible and reflect a superposition of almost filled and 
almost empty levels. There are no LZ transitions in that 
limit, but only a continuous change of energy bands at 



the branching times where E = E ^. However, for 

larger V/T, the LZ probability becomes sizeable and the 
dynamics is more complex, generally involving a dynam- 
ical population of all four subgap states. The right panel 
in Fig. [2] displays the case of relatively large V, where the 
system oscillates due to LZ transitions between the levels 

and E^ The frequency ojs of the S{t) oscillations 

is much slower than the Josephson frequency loj = 2eV/h 
and determined by the lowest interlevel transition energy, 

LUs = min(£^ |_ — E ). Note that cos reappears in the 

ac Josephson current. 

To conclude, we predict that the conformational degree 
of freedom (represented by a TLS) in a superconducting 
molecular dot or break junction responds in a dissipa- 
tionless manner to variations of the phase difference ip 
across the dot /junction, including a complete reversal. 
This effect should be observable using existing experi- 
mental methods over a wide parameter range. Under an 
applied voltage, this effect leads to quasi-periodic TLS 
dynamics due to the time-dependent occupation prob- 
abilities of Andreev states. — We thank T. Martin for 
discussions. This work was supported by the SFB TR 
12 of the DFG and by the EU networks INSTANS and 
HYSWITCH. 
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